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I. INTRODUCTION 

Now that the a'^{Za) correction to the hyperfine sphtting of the ground state muonium is 
known very accurately because of recent works ||l|-3, the previously-calculated a(Za)^ term 
has become one of the main sources of the remaining theoretical uncertainty. Improvement of 
this error is urgently needed in view of the new muonium hyperfine measurement in progress 
Q. Here Z is the "atomic" number of the nucleon. Of course Z = 1 for the muonium, but it 
is kept to indicate the bound state origin of the correction terms. In this paper, we present 
the a(Za)^ radiative correction evaluated in the NRQED bound state formulation |^. 

The previous evaluation of this term, including recent modifications H,^ , gives 

a{Zay 



Au[a{Zay] = Ep- 



TC 



3 ^ ' 
37 8 



|ln4 + ^ + ^)ln(Z«)-i + 14.88(29) 



(1) 



where the Fermi frequency Ep is defined by 

(2) 



^ 16 2 n "^ 

Ep = — a cKoo^rT 
3 M 
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Roo is the Rydberg constant for infinite nuclear mass, and m and M are the electron and 
muon masses, respectively. The presence of ln{Za) factors is due to the near infrared (IR) 
structure of the loosely bound state. The coefficients of In^ (Za)~^ and \n{Za)^^ terms 



were obtained by Layzer [g and Zwanziger independently. Brodsky and Erickson |[T0| 
confirmed these logarithmic terms and gave the leading contribution of the non-logarithmic 
term. Sapirstein reported the numerical evaluation of the non-logarithmic constant due to 
the radiative photon [|lT|. For convenience's sake, we refer to this non- logarithmic constant 
of the a{Za)'^ correction as the BES term. 

To compute the BES term, Sapirstein started from the relativistic bound state formalism 
and evaluated the entire a{Zay term numerically. In his approach, only the double loga- 
rithmic term was confirmed by varying the "atomic" number Z. Since the logarithmic term 
is a consequence of the near IR singularity of the bound state, the convergence of numerical 



integration worsens in the region of small momentum. The uncertainty in the BES term 
comes mainly from this difficulty in the numerical integration. Note also that his result 
contains terms of higher orders in Za. 

Our calculation of the BES term starts from the NRQED formalism proposed by Caswell 
and Lepage |]^. This approach enables us to isolate the a{Zaf term without being tangled 
up with higher order terms in Za: All these terms arise from different parts of the NRQED 
Hamiltonian. The leading logarithmic contribution is analytically separated. The small 
photon mass A is used in our approach but the A dependence can be easily identified and 
analytically subtracted in the numerical evaluation of each diagram. This is important for 
reducing the computational error of the BES term. 

Asides from the factor [aiZa)'^ j-n^Ep our results for the vacuum polarization contribu- 
tion and the radiative photon contribution to the BES term are 

A^BESvp = -0.218 567 ■ ■ ■ , (3) 

and 

A^BESph = 17.122 7 (11) , (4) 

respectively. Our final result for the total aiZaf correction is 

aiZa)' 



^v\a{Zay\ = Ef- 



n 



--\n\Za)-' 



+ ( -^ ln4 + ^ + A ) ln(Z«)-i + 16.904 2 (11) 



(5) 



The BES term has also been evaluated recently by Pachucki [|13| using the method he 
developed for the a{Za)^m Lamb shift correction |jl^. His result for the radiative photon 
contribution is 

AEBESph = 17.122. (6) 

Although this is in good agreement with our result (^, there remains some disagreement in 
the details. This will be discussed in Sec. VI. 



In Sec. II, we briefly describe our approach, namely the NRQED method, to the bound 
state problem. A more detailed prescription of NRQED is found in Ref. [Q. In Sec. Ill, 
the well-known Breit (Za)"^ correction ||T^ is rederived from NRQED. This calculation. 



multiplied by an appropriate "renormalization" factor, actually provides a part of the alZa)"^ 
term. It also serves as a prototype of the calculation of the entire a{Zay correction. In 
Sec. IV we derive the a(Za)^ contribution arising /^from the vacuum polarization insertion. 



We uncovered a mistake in the previous calculation in Ref. fl^ as was mentioned in our 
previous paper P]. The detail of the a{Za)'^ correction due to the radiative photon on the 
electron line is described in Sec. V. Sec. VI is devoted to the discussion of our result. 

II. OUTLINE OF THE NRQED METHOD 

The NRQED is a theory with a finite UV cut-off, which is completely equivalent to QED 
when it is applied to low-energy systems with typical momenta less than the UV cut-off A. 
The NRQED Lagrangian consists of operators which satisfy the same symmetries as the QED 
operators except that they satisfy Galileian invariance, although the final observable results 
of calculation are Lorentz invariant. Fermions are represented not by Dirac spinors but by 
Pauli spinors. The NRQED Lagrangian can be divided into two parts: The main part //main 
consisting of fermion bilinear operators multiplied by up to two photon operators or pure 
photon terms and the contact interaction part Lcontact involving four or more fermions. Both 
parts of the Lagrangian are determined by the following simple rule: The operators which 
appear in its Lagrangian and their coefficients are chosen so that any scattering amplitude 
calculated in the NRQED coincides with the corresponding scattering amplitude of the 
original QED at some given momentum scale, e.g., at the threshold of the external on-shell 
particles. This matching condition is applied order by order to the expansion in the coupling 
constant a and velocity v of the external fermion. The Coulomb gauge is used in the NRQED, 
while the Feynman gauge is more convenient to compute the QED scattering diagrams. 
Readers interested in NRQED may refer to Refs. PJT^|-P^ . The precise description of the 



NRQED Hamiltonian can be found in Ref. [Q. After determining all operators and their 
coefficients to the desired order of velocity v of the electron and the coupling constant a, we 
evaluate the energy shift, etc. using the Rayleigh-Schrodinger perturbation theory, choosing 
as the unperturbed system the exact solution of the nonrelativistic Schrodinger Coulomb 
system. 

The main part of the NRQED Hamiltonian -f^^ain needed to compute the a^ZaY cor- 
rection terms is of the form 



H„ 
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where p and p' are incoming and outgoing electron momenta, respectively, and g'is incoming 
photon momentum. ^ The superscript A in -ffmain indicates that the theory is regularized by 
the UV cut-off A for the radiative photon. A may be less than m. The "renormalization" 
coefficients in (0) are 



2ap, 
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^ The term qi x A{qi) in Eqs. (53) and (55) of Ref. Q must be replaced by q2 x A{qi 



Cp = — . (8) 

The first term in (^ is the nonrelativistic kinetic energy term. The rest are named succes- 
sively as Coulomb, p^ relativistic kinetic energy, dipole coupling, Fermi, Darwin , seagull, 
iy-(wave function) derivative Fermi, q^ derivative Fermi, and p'p coupling, respectively. 
The last two terms bilinear in photon operators are introduced to represent the vacuum 
polarization insertion in the transverse and Coulomb photon propagator, respectively. 

For the muon line, only the Coulomb and Fermi terms are needed for the calculation of 
the a{ZaY correction. They are obtained by replacing e and m of the electron interaction 
terms by —Ze and M, respectively | . 

As for the contact part of the NRQED Hamiltonian -ff intact 5 only the spin-flip type is 
needed: 

where x is the Pauli spinor for the positive muon and di can be written as 
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di = a(Za)^7r- 

+ a{Zaf'^ 



■ ^ 13 3 
In 2 h - 

4 4 



rf(3) J— ^ + d(^) In^ f A) + rf(i) In f A) + rf(0) 
y [[Za)m)'^ \mj \mj 



[Za)m \m 



(10) 



Here A' is the UV cutoff of the transverse exchanged photon momentum. Note that, when 
the a{Zaf' part of Eq. ([T0|) is evaluated in the first order perturbation theory, it gives the 
a{Za) correction to the hyperfine splitting of the ground state muonium calculated by Kroll 
and Pollock, and Karplus, Klein, and J. Schwinger [^. For brevity, let us call the a{ZQ) 
correction as the KP correction. The coefficients d^^\ ■ ■ ■ , d^^"* of the a{ZaY part of Eq. (p!OD 
are pure numbers. These are the quantities that we want to calculate in this paper. 



^We use the convention that the electron charge is e and the positive muon charge is —Ze. 



The Green function Go{p, q', E) appearing in this calculation is known in an exact closed 
form for the nonrelativistic Coulomb potential [^ . For an arbitrary energy E in the complex 



E plane this function takes the form: 

2m 



Go(p, g; E) = -^^-^{27ry6\p- q) 
—2m —AirZa —2m 



GAtt ~ 

R{p,q;E) , (11) 



^^{Za)' 



where 



and 



R{p,q;E) 



pI = -2mE, u=^, (12) 

Po 
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[p'^ + Po){q^ + Po) 



"" Jo "^^Mp- qlVo + (1 - Py{p+Pl){q'+Pl) ■ ^^^^ 



The first, second, and third terms of (|TTD correspond to zero, one, and two or more Coulomb- 
photon exchanges. For E = Eq — k, where Eq is the ground state energy and k is the energy 
of a radiative photon, pg may be written as 

pI = -2m{Eo - fc) = 7^ + 2mk . (14) 

We calculate the {Za)^ and a{Za)'^ corrections in the subsequent sections using the NRQED 
Hamiltonian and the Green function given above. 

III. THE {Zaf CORRECTION 



In this section, we rederive the well-known Breit {Za)^ relativistic correction |T^ from 
the NRQED in order to illustrate how it works, particularly, how the contact term ifcontact is 
constructed. The first computation of the Breit term in the framework of the NRQED was 
carried out in Ref. f^. In the previous paper [Q, we have shown that both a{Za) and a^{Za) 



corrections come from the first order perturbation theory of the -ff^ntact which represents the 
difference between the QED scattering amphtude and the NRQED scattering amphtude. In 
contrast, the contact term calculated here is derived entirely from the NRQED scattering 
amplitudes. Their presence is crucial to cancel the UV divergence occurring in the bound 
state calculation of the operators in the -ffmain- More importantly, the calculation carried 
out in this section immediately yields a part of the a{Za)'^ correction if one takes account 
of appropriate "renormalization" constants of the potentials. We must make sure that this 
calculation is consistent with the calculation of the other parts of the a(Za)^ correction. 
This is why we want to rederive the Breit correction by the NRQED method, not by other 
bound state formalism. 

A. Diagram Selection 

The first task is to identify the potentials contributing to the Breit (Za)"^ correction. 
This can be achieved using the order of magnitude estimate of various operators appearing 
inthei/li^ofEq. ® 



< d >r^ m{v/c), < 9t >~ m(i;/c) , < eA" >~ m(t;/c) , 

< eA >~ m{v/cf, < eE >~ m'^{v/c)^, < eB >~ m^{v/c)^, (15) 

where (f/c) ~ (Za). For example, for the Coulomb potential between the electron and the 
muon 

^c(p',p)--^, (16) 

where k = p' — p, we have < Vc > /m ~ {v/cY ~ {Zaf. We will set the photon mass A to 
zero in the bound state calculation. 
Since the Fermi potential 

, -iekxa^ iZe{-k) x a -1 



has an expectation value of order {Za)'^{m/M)m, one source of the Breit (Za)^Ep correction 
is the first order perturbation with the order {Zaf'{mjM)m potentials. All potentials of 
this type must have spin flipping property in order to contribute to the hyperfine structure. 
One of these potentials is the P^-derivative Fermi term in Eq. (^ which yields the potential 



Vw{p\p} 



8m3 2M fc2_^A2 



Another contribution comes from the seagull term: 

, ^ ie^qxcfe iZe{-k) x a^ "1 / v \ ^ 



4m2 2M A: 2 + A2 g ^ + A^ ' 

— * 

where k and q are the transverse and Coulomb photon momenta, respectively. 

One must also consider possible contributions from higher order perturbation theory. 
The second order perturbation term has the form 

E - i^„=i 

_ V- {'^n=l\V\'^k){'ipk\V\'^n=l) .^„x 

Ml ^n=l ^k 

Since the denominator E^ — E^ is of order {Za)'^'m, one potential must be of order 
{Zay{m/M)m while the other is of order {ZaYm in order that 6E contributes to the 
Breit correction: 

6E ~ iZ<^rMM)miZarm ^ ^^^^a - ^ _ (21) 

This can be realized only if one is the Fermi potential and the other is an order (Za)^m 
spin-non-fiip potential. We find two candidates for the latter: the relativistic kinetic energy 
term 

\3r3/ 



VK{p',p) = VKm'^'nY6\p-p' 

\3r3 



^^{27.f5\p-p') 



im-^ 



1 
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{Vk) ~ ^ ~ {Zafm , (22) 



rn^ 



and the Darwin term 



Vn{p',p)^^{-Ze)- ^ 



m2 ^ 'A;2 + A2 
(Vn) ~ ^^^ ~ {Za)'m . (23) 

The third and higher order perturbation terms do not contribute to the Breit correction. 

B. Determination of the NRQED Contact Terms 

We have identified the NRQED potentials necessary for the calculation of the Breit cor- 
rection, namely, the Fermi, derivative Fermi, seagull, relativistic kinetic energy, and Darwin 
terms. The next step is the determination of NRQED contact terms corresponding to these 
potentials. The contact terms of spin-non-fiip type contribute to the hyperfine splitting 
calculation only through the second- or higher-order perturbation terms, analogous to the 
case of the Darwin or relativistic kinetic energy potentials. However, there is no {ZaYm 
contact term of the spin-non-fiip type. The lowest order spin-non-fiip contact term is of 
order a{ZaYm, which gives the relativistic binding correction to the Lamb shift. Therefore, 
we have only to consider the NRQED contact terms of spin-fiip type 

SH = -d,-^ii,^a,ij) ■ ix^a^x)- (24) 

In the following the muon is treated as an external static photon source since the infinite 
muon mass limit is taken. However, the method of the contact term determination for the 
dynamical muon case is essentially the same as the static case and the latter can be readily 
adapted to the former. 

As was explained in Ref. |^, the contact terms are determined by comparison of QED 
scattering amplitudes and NRQED scattering amplitudes. The lowest order diagrams con- 
tributing to the contact term relevant to the Breit correction is the two-photon-exchange 
process between the electron and muon, namely, one-loop process. In order to contribute 
to the Breit correction, at least one of the photons must be transverse. Exchange of two 

10 



transverse photons results in a recoil type correction proportional to {m/MY and is not of 
interest here. That leaves diagrams with one transverse photon and one Coulomb photon. In 
QED, there is one diagram with 67* and 67° vertices. NRQED interaction terms which give, 
in the hound state calculation, higher-order contributions than the order we are interested in 
are to be ignored in the scattering amplitude comparison. This reduces the relevant NRQED 
scattering amplitudes to the following five combinations: a Coulomb potential with a Fermi 
potential, a Fermi potential with a relativistic kinetic energy term and a Coulomb potential, 
a Fermi potential with a Darwin potential, a Coulomb potential with a derivative Fermi 
potential, and a seagull potential. The first four are given by the second order perturbation 
theory of the NRQED Hamiltonian and the fifth is from the tree NRQED Hamiltonian. 
These five together determine the contact terms represented by the shaded circle in Fig. 
|l](a). Only the p^ /2m part of the NRQED Hamiltonian is treated as unperturbed system 
for the scattering perturbation theory. In other words, the Coulomb potential appears as 
one of perturbative potentials in this comparison. The comparison are shown in Fig. |l](a). 

In this case the QED scattering amplitude at threshold is completely replicated by the 
NRQED scattering amplitude consisting of a Coulomb potential and a Fermi potential. Thus 
the contact term must be chosen as minus sign times the sum of the remaining four NRQED 
scattering amplitudes. 

The one-loop contact terms determined in this comparison are all linearly divergent 
and hence their values depend on how they are regularized. Although gauge invariant 
regularization method is desirable, it is possible in this case to use a simple momentum cut- 
off. This is because, as we shall see in Appendix A, calculation of the hound state expectation 
value also leads to a divergent integral and must be regularized. Even though the regularized 
bound state calculation is not gauge invariant, the gauge invariance can be restored if one 
chooses appropriate contact terms. What is crucial is that the regularization method of the 



contact term is consistent with the bound state calculation ||22[ . We satisfy this requirement 
by introducing an UV cut-off A in the momentum of the transverse exchanged photon in 
both contact term calculation and bound state calculation. 

11 



The NRQED scattering diagrams can be easily written down using the NRQED Feynman 
rule given in Fig. 3 of Ref. [^. Against the relativistic kinetic energy term this rule leads 
to the contact term of the form 

^le-. = -2 t^V,{Q, k)^VK{k)^Vc{k, 0) . (25) 

Here the factor two is to take account of the time-reversed diagram. The contact terms 
against the Darwin term, derivative Fermi term, and seagull term can be similarly con- 
structed and evaluated. Explicit evaluation of (^) and other terms is carried out in Ap- 
pendix A. 

Next we consider the three photon exchange process, or two-loop one, shown in Fig. 
|I|(b). This requires adding one more Coulomb photon exchange potential to both QED and 
NRQED amplitudes of the two photon exchange type. Adding another type of potential of 
NRQED, such as the Darwin potential, is not necessary since it gives a power of v"^ ~ (^«)^ 
higher than what we are interested in. However, it is necessary to consider new kinds of 
scattering amplitudes, which are the combinations of the contact terms introduced in the 
two photon exchange comparison and one Coulomb photon potential. Again, the QED 
(left-hand-side of Fig. |I|(b)) contribution is equivalent to the first diagram of the NRQED 
contribution (on the right-hand-side of Fig. |l](b)). Thus the contact terms are determined by 
the rest of the NRQED scattering amplitudes. Processes with four or more photon exchange 
clearly have more explicit powers of Za, yielding higher order contact terms. Thus we don't 
have to deal with them as far as we are interested only in the (Za)"^ Breit correction. 

Some of the two-loop contact terms, which have a Coulomb photon exchange potential 
sandwiched between two other kinds of potentials, are logarithmically divergent in both UV 
and IR regions and will kill UV divergences in the bound state calculation. Others are 
linearly divergent in both UV and IR. These diagrams have the structure such that one 
Coulomb photon potential is added to the edge of one loop linear UV divergent diagrams 
proportional to A/m introduced in the two photon exchange comparison. The additional 
Coulomb photon on the edge yields the threshold IR singularity causing the term to be 
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proportional to m/A, where A is the photon mass introduced as the IR cut-ofF. The multi- 
phcation of the UV and IR divergences results in the form A/A. The UV-IR singularity A/A 
appearing in these diagrams is completely canceled out by the diagram which consists of 
the one-loop contact term introduced before and the Coulomb photon potential connected 
by the free fermion propagators. 

Here we show one of the two-loop contact terms against the relativistic kinetic energy 
term: 

Other contact terms related to the relativistic kinetic term in this order are presented in the 
Appendix A. The potential V2c-k(i) is logarithmically divergent. Other terms are linearly 
divergent. 

The contact term against the combination of the one-loop contact term and the Coulomb 
potential, which we denote with the suffix lloop^, is also linearly divergent and is given by 

V2c-Kiioo,^ = ^ j j^, J j^,Vp{0,k)^^VK{k)^^Vcik,0)^Vc{l,0) . (27) 

Other two-loop contact terms involving the Darwin, derivative Fermi, and seagull terms are 
given in Appendix A. 

C. Summary of the (Za)"^ Correction 

We have prepared the non-radiative NRQED Hamiltonian, including the contact terms, 
up to the order {Za)^{m/M)m: 

-"NRQED — r 

eAf)^. (T„ ■ m p^n- Ax E - 



2m r 2m 8m^ 8m^ 

e{p2, (fe- B} e^a ■ Ax E 



]m^ Am? 



- ^i;;^(V^^^"^eV^) ■ (^^^7.^) • (28) 

The contact term coefficient di is calculated in Appendix A. It is the sum of the contributions 
from the relativistic kinetic energy, Darwin, derivative Fermi, and seagull terms: 
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di = {Zaf-Ti 



-2 21nf-) -41n2 + 61n3 

+— + lnf-) +21n2-31n3 + - 

77r V A / 2 

+2 — + 41n2-41n3 

77r 

A , /A\ , , 1 
+ ln T -21n2 + ln3-- 

77r VA/ 2 

. 



(29) 



The value of di could vary for different regularization methods. Although di adds up to 
zero in our regularization method, this does not mean that we do not need the contact term. 
Finiteness and gauge invariance of the final answer are guaranteed by the presence of this 
contact term in individual terms. 

Using the non-relativistic Coulomb system as the unperturbed system, we can now calcu- 
late the binding effect, namely, the Breit hyperfine energy correction in perturbation theory. 
The results are summarized here as terms proportional to the Fermi energy Ep. AE^, AE^, 
AEyj, and AEg are the contributions from the relativistic kinetic energy, Darwin, derivative 
Fermi, and seagull terms, respectively (see Figs. ^ ^ ^ and|^): 

1 



AEk = Ep^Zaf 

AEd = Ep^Zaf 
AE^ = Ep{ZaY 
AEs = Ep^Zaf 



--61n(f)+21u ^ 



3.„|1 



< 



- + 4 In f - 



(30) 
(31) 

(32) 
(33) 



where 7 = Zam is the typical momentum scale of the muonium in the infinite muon mass 
limit, and A is the photon mass which is set to zero at the end of the calculation. The sum 
of Eqs. dg), (0), (|2[), and (§D gives 



Av[{Zaf] = -EpiZa) 



(34) 
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which is the well-known Breit correction. Note that the IR singularities ln(A/7) cancel out 
when all diagrams are summed up. In Appendix A, we show details of derivation of these 
terms. 

We have shown how the NRQED bound state formalism works using the well-known Breit 
correction to the muonium hyperfine structure as an example. As we have seen, the contact 
term in NRQED plays the crucial role: it describes the high energy behavior, recovers the 
symmetry, such as Lorentz symmetry and gauge invariance, and kills the would-be IR- and 
UV-divergent quantities. 

IV. THE a{Zaf VACUUM POLARIZATION CONTRIBUTION 

The order a radiative correction in the a(Za)^ term comes from two sources: One is from 
the vacuum polarization insertion in one of the exchanged photons between the electron and 
muon and the other is from the spanning photon on the electron line. They are separately 
gauge invariant. In this section, we will deal with the contribution coming from the vacuum 
polarization insertion. 

The result of our numerical evaluation of the vacuum-polarization contribution was 

a{Zaf 



Aiy[a{Za) 



Jvp 



-Ef 



n 



Inf— j -0.218 81 (29) 



(35) 



15 \7 

where the error comes from the numerical integration. The error associated with the finite 
photon mass, which was used as an IR regulator, is of order A/m, hence negligible for the case 
X/m = 10^^. Our result Eq. (^) disagreed with that obtained by Brodsky and Erickson 

in 
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-0.085 234- ■• 





(36) 



This term actually consists of two parts. One corresponds to the vacuum polarization 
insertion in the Coulomb photon, and the other is insertion in the transverse photon. Our 
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result for the first part agrees with the corresponding result of Brodsky and Erickson. For 
the second part, however, we found —0.801 (4) instead of —2/3 in Eq. (pGl). In an effort to 
determine the cause of the discrepancy, we have analytically evaluated the integral expressing 



the vacuum polarization contribution, given by Zwanziger in Ref. |2^. Our subsequent 
analytic work showed that the contribution for insertion in the transverse photon is —4/5, 
in agreement with our numerical result. Using this corrected value, the numerical value of 
the non-logarithmic part of Az/[a(Za)^]vp becomes 



aiZoi 



i2 



AE'besvp — Ep 

7r 



8 , 214 4 

-15 ^"' + 225-5 



"^^"^-^^(-0.218 567 ■■■) . (37) 



TT 



Since we found the error in the calculation of Ref. JTOl (in June 1994), Sapirstein has also 
noticed it independently ||2^, and Brodsky and Erickson have agreed with the corrected value 



| 25| . We were also informed by Karshenboim |2^ that the same result was independently 



obtained by Schneider, Greiner and Soff |^ 



Let us briefly describe the NRQED evaluation of Eq. (|35|) . The terms relevant to the cal- 
culation of the NRQED Hamiltonian (0) are the Coulomb, Fermi and two-photon interaction 
terms. Combining these three interaction terms, we can construct the effective spin-flip and 
spin-non-fiip potentials: 

and 

Kvp - -^^p(V^V)(x^x)(^4^ ' (^^) 

respectively. We will construct the contact term of the NRQED Hamiltonian of Eq.(^ 
by comparing the QED and NRQED scattering amplitudes using the same recipe as for 
the Breit correction discussed in Sec. II. In this case both QED and NRQED scattering 
amplitudes contribute to the contact term. (See Fig. ^). We evaluate the effect of these 
potentials using the nonrelativistic Rayleigh-Schrodinger perturbation theory. 
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To see the UV and IR cancellation of the calculated result, however, it is more conve- 
nient to decompose the contribution of -ff^ntact i^^^ several parts. The detail is described 
in Appendix Bl. In this treatment, the bound state contribution comes from the first-order 
perturbation of the spin-flip potential Vi;vp of Eq. (^) and from the second-order pertur- 
bation involving the spin-non-flip potential Kvp of Eq. ( PD| ) and the Fermi potential Vp. 
The UV divergences in both bound state calculations are taken care of by the corresponding 
contact terms derived from the scattering amplitudes of these NRQED potentials. The sum 
of two contributions is found to be 



2 



A 771 a{Za) „ 

^-C/NRQEDvp = ^Fzr: 

71 15 



In(^')_91n2 + -ln3-- 
7/ 2 4 



(40) 



The remainder is a part of the -ff^ntact coming from the QED scattering amplitude. The 
three-photon-exchange diagram with one vacuum-polarization insertion contributes to this 
order. Since, at the two-photon-exchange level, we have introduced the contact term in 
NRQED Hamiltonian representing the vacuum polarization effect of the a{Za) KP term, 
which is the — 3/4 term on the first line of Eq. (|TU|), we must also take into account the 
contribution of this term in evaluating the NRQED scattering amplitude. The NRQED 
scattering diagram with the Kroll-Pollock (KP) contact term connected with the Coulomb 
potential by the free electron propagator should be subtracted liiom the QED scattering 
amplitudes. Their numerical value, for the photon mass X/m = 10~^, is 



A£;qedvp = -^—LEf{5.520 74 (29)) 
aiZal 

= rL/p 



n 



-^ln(^) -0.619 49 (29) 
15 V A 



(41) 



The sum of AE'nrqedvp and A_Eqedvp gives the result of Eq. (^). The cancellation of the 
photon mass dependence between AE'nrqedvp and AE'qedvp has been checked numerically 
using various values of A. Of course Eq. ( P^D is superseded by the analytic result which 
incorporates (P7|). 
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V. THE a{Zaf RADIATIVE PHOTON CONTRIBUTION 

More challenging is the radiative photon contribution to the a{Za)'^ term. The sources 
of the a(Za)^ correction have been discussed in Ref. 0. Starting with the main part of the 
NRQED Hamiltonian (|^), we can construct the contact term -ff^ntact- Then we find four 
types of contributions to the a(Za)^ correction in the bound state perturbation theory. The 
manifest cancellation of UV and IR divergences occurs between the bound state calculation 
of the operators from the -f^main ^^^ ^^^ corresponding part of the contact term derived 
from the same operators. For simplicity we omit the overall factor {a{Za)'^/TT)EF in the 
following. This means also that the anomalous magnetic moment Oe stands for 1/2. 

One contribution, denoted Eb, where B implies the bound state effect, is the third order 
perturbation with two dipole couplings and one Fermi potential from the main part of the 
NRQED Hamiltonian, -f^main- (See Fig. ^) The radiative correction of order a comes from 
the presence of an intermediate virtual photon. For this calculation, we use the Coulomb 
Green's function in the presence of a photon with energy k obtained /^from Eq. ([T3|). The UV 
cut-off A is imposed on the radiative transverse photon momentum. All UV divergent parts 
and some finite parts are analytically calculated. Remaining finite parts are numerically 
evaluated using VEGAS [^]. The results are listed in Table |. All relevant formulas are 
given in Appendix B2. 

The second contribution Eg comes from the contact interaction part of -^contact corre- 
sponding to the bound state calculation of Eb- The NRQED scattering diagrams with two 
dipole couplings and one Fermi potential yield this contact interaction term in NRQED. (See 
Fig. 1^.) The UV cut-off A is also imposed on the radiative transverse photon momentum. 
The calculation is entirely analytic. The results are listed in Table |I|. 

The third type of contribution E^,, where R stands for the "R" enormalization constant, 
is found in the first and second order perturbation terms with the potentials which give 
the {Zay correction. These potentials are found using the analog of the method of the 
Breit (Za)'^ correction. In this case, the radiative correction of order a comes from the 



"renorinalization" coefficients of order a of these potentials. The UV cut-off A' is imposed 
on the exchanged transverse photon momentum. This A' cancels out when the bound state 
calculation is completed. Some of these contributions are already calculated in Sec. Ill 
concerning the Breit correction. For the contributions from the relativistic kinetic energy, 
Darwin, and seagull term, we only need to multiply Eqs. (pOD, (PT|), and ( PB| ) with the 
appropriate "renormalization" constants given in Eq. (||). The diagrams for these contri- 
butions are shown Fig. H, ^,and ^, respectively, and the formulas are given in Appendix A. 
For instance, the contribution associated with the Darwin term is given by 

A^rf = (cb cf - l)AEd, (42) 

since it involves both Darwin and Fermi term. The remainder, the contributions associated 
with the g^ derivative Fermi and p'p interaction term, requires new calculation. The diagrams 
of the g^ derivative Fermi term are analogous to the derivative Fermi term shown in Fig. ^. 
The diagrams of p'p interaction term is shown in Fig. ^. The formulas for both contributions 
are listed in Appendix B2. The result of Eji evaluation is summarized in Table |T|. All terms 
proportional to In (Za) and ln(Za;) can be obtained from E^ and Eji. We find 



,2- 



which agrees with the known result by Layzer |Q and Zwanziger p. Here, the ln(?7i^/7^) 
term is related to the IR behavior of the radiative photon and the ln{m/'y) term comes lirom 
the threshold singularity due to the exchanged Coulomb photon. The double logarithm 
In (m/7) is a consequence of simultaneous presence of both types of contributions. 

The last type of contribution Eq comes from the contact term of the NRQED Hamil- 
tonian which is calculated from the QED scattering amplitudes exchanging three photons 
between the electron and the muon, dressed by one radiative photon on the electron line. 
/^From this one must subtract diagrams which are the KP contact term combined with one 
Coulomb potential and the Fermi potential times anomalous magnetic moment combined 
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with two Coulomb potentials. (See Fig. ^) The numerical results evaluated by VEGAS 
for various photon mass are listed in Table |V| . The IR regulator mass A is kept finite in the 
exchanged photon in the scattering diagram calculation. The radiative photon is chosen to 
be massless. This is allowed since the radiative and exchanged photons are separately gauge 
invariant. The QED scattering contribution depends only on m/\ and may be written as 

Eq = ci In^ (^j^ + C2 ln(^y) + C3 + C4- + ■ ■ • . (44) 

The coefficients ci, C2, ■ ■ ■ can be determined by fitting the formula with the result of 
numerical integration pS]. Ignoring C4, C5, ■ ■ ■, we found 



ci = -2.664 4 (4), C2 = 1.171 5 (71), C3 = -8.530 (28). (45) 

Actually, the logarithmic part 

(EQ)iog = Ciln2(^^)+C2ln(^) (46) 

can be determined analytically, noting that the sum of the contributions Eb, Es, E^i, and 
Eq must be free of both UV cut-off A and IR cut-off A. We find that 

{EQ)iog = -{Es + Epilog 

^From this we obtain 

Ci = — = —2.666 666 • • • , 
3 

11 7 

C2 = 241n3-321n2 Oe = 1.213 762 9 ■■■ , (48) 

9 2 



which are in good agreement with the numerical result (^5D.p| Using the analytic result (^8]) 
for ci and C2, we can determine the constants C3 and C4 more precisely from the numerical 
result. We find m 



Hi we choose a = -8/3 in Eq. (H), we obtain C2 = 1.208 9 (5) and C3 = -8.675 8 (46). 
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C3 = -8.724 3 ± 0.001 1, C4 = 8.688 ± 0.825 . (49) 

Using the constant C3 in Eq. ( ^9|) together with other contributions summarized in Tables 
I, II, and III, we arrive at the final result for the radiative photon contribution to the BES 
term: 

A^BESph = 17.122 7 (11) . (50) 

VI. DISCUSSION 

The value (6) of the a(Za)^ correction obtained by Pachucki [13] is in good agreement 
with our result (4). However, it is not fully satisfactory in the sense that details of his 
calculation seem to disagree with ours beyond the uncertainty of numerical evaluation. The 
numerical value of his low energy contribution Ei^ (defined by Eq. (16) of [13]) is smaller 
than our E^ by 0.025 7 (2). Unfortunately, his mid energy contribution Eu or high energy 
contribution En cannot be compared directly with our Eg-, E^, or Eq. We find, however, 
that the sum Em + Efj is larger than our sum Es + Er + Eq by 0.024 6(14). Note that 
these discrepancies nearly cancel out, leading to a good agreement of (4) and (6). 

In order to compare these results in detail, let us examine El more closely. E^ is 
the low energy contribution characterized by the parameter e which is the (UV) cut-off of 
the frequency u of the virtual photon. In Ref. [13], terms of E^ divergent for e — * 00 are 
separated out and evaluated analytically. The remainder is evaluated partly analytically and 
partly numerically. Adding up the results (24)- (27) of Ref. [13] gives the sum of these terms 
^from which the contribution of mass renormalization term is not subtracted. Meanwhile, 
the mass renormalization effect is subtracted out in Eqs. (B19), (B20), and (B21). 

One way to compare these two methods is to drop the mass renormalization terms from 
(B19), (B20), and (B21). Then evaluate them, for instance, for an NRQED UV cut-off A 
satisfying Am/'-f'^ = 1. Correspondingly, we choose em/'y'^ = 1 in Pachucki's formula for El. 
We have evaluated both of them numerically. The agreement between the two integrated 
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results is excellent: NRQED gives 0.455 763 (0.000 033), while Pachucki's formula gives 
0.455 769 (0.000 010). 

In order to identify the leading contributions to the alZa)"^ correction we have also 
performed an asymptotic expansion of our formulas (without mass renormalization terms) 
in 7^/(A?7i). The corresponding procedure for Pachucki's E^ is to take the limit of small 
7^ first and then go to vanishing e. Both procedures are found to give the same result for 
the logarithmically divergent and finite constant terms. This is as expected since the mass 
subtraction term in NRQED has a linear divergence proportional to Am/'-f'^ without any 
constant term. Q Note that the corresponding term in Pachucki's method is proportional 
to em/'-f'^. Thus the resultant logarithmic and constant terms must be of the same form in 
both Eb and E^, if we identify A = e. These checks have convinced us that El and Eb are 
exactly equivalent as far as their analytic properties are concerned. 

As the second check, we have also evaluated Pachucki's formula for El, after carrying 
out an explicit mass renormalization, by numerical integration for values of A(= e) ranging 
from 10^'^7^/?7i to 10^'y'^/m and extrapolating the result to A = oo. This method is adopted 
because precision of numerical integration deteriorates steadily with increasing A and direct 
evaluation of our formula at A = oo becomes impossible. In spite of our considerable 
computational effort over several months, however, the resulting El — Eb obtained thus far 
is not capable of distinguishing between the values and the difference 0.025 7 mentioned 
at the beginning of this section. 

In conclusion, we have seen nothing wrong thus far with the analytic form of either El 
or Eb but have been unable to resolve the apparent numerical discrepancy. Neither can we 
rule out an intriguing possibility that Pachucki's separation into low and higher energy parts 



^In fact, in the NRQED approach there is also a contribution proportional to ^ tvaij^^ from the 
diagram i?oo of Fig. 7. This term is cancelled out by the same but negative contribution from the 
diagram 5o of Fig. 8. 
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is slightly different from ours but the sum is not affected by how the separation is made. 

We publish our result in the present form to avoid undue delay and to encourage further 
investigation by others. 



The result obtained by Sapirstein for the radiative photon contribution was 11 



A^BESph = 15.10 (29). (51) 

Concerning the apparent discrepancy between Eqs. (^ and (|5T|) , Karshenboim pointed out 
that the latter might include the known a{Za)^ln[Za)^^ correction [^,^ 



EpaiZaflnf-) (5\n2-j-) = E^a(Za)Mnf-V-8.471 764 ■ ■ ■) 

= Ef— ^(-0.955 ■■■) 52) 

TT 



in the relativistic bound state formalism adapted in Ref. [TTI]. We wrote the second line 
of (^) to emphasize that, numerically speaking, this contribution is of the same order of 
magnitude as the a{Za)^ term. 

Eq. (^) will contain in addition the non-logarithmic part of the a{Za)^ correction. If 
terms of order a{ZaY or higher can be ignored, this contribution may be estimated from 
Eqs. (§, (|1]), and (0) to be 

EpaiZaf {-37 ±13). (53) 

In order to see whether this is plausible or not, let us note, based on the NRQED analysis, 
that the non-logarithmic term of the a{Za)^ term may be smaller by only a factor of 
In I — j ~ 5 compared with the leading term Eq. (|52|) : 

/ 191\ 
EpaiZaf i5\n2 j x C(l) = Epa{Zaf{-SA71 764- ■■) x 0(1). (54) 

Thus the estimate (|53|) is somewhat large but not inconsistent with our NRQED estimate. 
The correction due to one radiative photon to the hyperfine structure was computed 
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recently for various atomic number Z p2| , |33| without expanding in powers of Za.0 Thus 
their results ioi Z = 1 must contain corrections of order a;(Za)" with n = 0,- ■ ■ , oo. Un- 
fortunately, the results at Z = 1 do not have enough accuracy in order to extract the 
non-logarithmic contribution of the a(Za)^ correction, namely the BES term. Their results 
for high Z are not sufficiently accurate either for a reliable extrapolation to Z = 1. 

In this circumstance it seems that only a direct evaluation of the pure a{Za)^ term will 
clarify the ambiguity in the theoretical calculation concerning this term. This problem will 



be treated in our subsequent paper [|J]. 

To summarize, we have calculated the a{Zay correction to the hyperfine splitting to 
the ground state muonium based on the NRQED. Including the logarithmic terms, the total 
contribution of this order becomes 



Ai'[a{Za) 



a{Za) 



-E, 



TC 



— In 

15 V7 



(t)- 



8 , o 214 4 

— In 2 H 

15 225 5 



--In^ 
3 



Ef 
+ 



m 

8 11 7 

1 — a, 

3 24 2 

ai^Za)"^ r 8 2/"^ 

vr [3 V 7 

8, . 37 8 

-ln4H \ 

3 36 15 



In 2 + 



Inf— j +17.122 7 (11) 



^^^{—^ +16.904 2 (11) 



(55) 



where Oe = 1/2. Note that the uncertainty due to the numerical integration is now reduced 
drastically to match the expected accuracy of the forthcoming experiment. 

Now that we have completed the evaluation of the a"(Za)'^~"', n=l,2,3, terms, it 
seems about time to discuss the numerical value of the theoretical prediction of the hy- 
perfine splitting of the muonium. However, the leading logarithmic corrections of order 
a^{ZaY~^,n = 0,1,2,3, turn out to be numerically of the same order of magnitude as 



^We thank P. Mohr and B. Taylor for drawing our attention to Ref. [^ and discussing their 
results. 
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some a^ corrections. Thus, we postpone comparison with experiment to the next paper in 
which we will treat these higher order logarithmic corrections. 
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APPENDIX A: CALCULATION OF THE {Zaf CORRECTION 

1. List of contact terms 



In addition to the contact term (25) against the relativistic kinetic energy term, there 
are also contact terms against the Darwin term, the derivative Fermi term and the seagull 
term determined from the two photon exchange scattering amplitudes. They are given by 

-^^VFiO,k)^^Vn{k,0), (Al) 

^Vw{0,k)^^Vc{k,0) , (A2) 

/A ^3/j, _^ 

j^,Vs{0,0,k), (A3) 
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where the potentials Vf,Vd,Vw,Vc, and Vs are defined by (|T3), (|2|), (0), (y), and (H), 
respectively. 

iFiom. the three photon exchange scattering amplitude, we find, in addition to (^), the 
following contact terms against the relativistic kinetic term: 

y^^-m = -2/ (^y (^^HO,A:)^^FK(fc)^^Vb(A:,0-^^c(/,0) , (A4) 

r^ d?k f (fl , r 7t -2m - - -2m 

V^2c-fc(3) = -2 / 77-^ / 7T^-T^V^c(0, k - l)-^^^VK{k - l)- 



2nrJ (27r)3 ^^ ' ^|^_/|2 "^ ^|^_;|. 



X \/^(A; - /, l)^^Vcil, 0) . (A5) 



2. The Relativistic Kinetic Energy Term 

The correction coming from the relativistic kinetic energy term shown in Fig. ^ is 
calculated in second-order perturbation theory given in (|2^). 

The diagram of Fig. ^(a) with two or more photons exchanged leads to a finite result 



AE[Fig.2(a)] 



2 



(27r)12 (^-2+ ^2)2(^2 +y) 

-* /— 647r\ ~ -* 
X {VF{f,l))[-^^,)R'{l,q;Eo)VK{q,p) 



TT JO (g^ + 7^) 
1 a^ 472 



ln2-- + -tan~^- 
2 q 7 



2 7^ g2 _)_ ^ 

= E,{Zaf^^ , (A6) 

where the factor 2 accounts for two contributing diagrams. The photon mass A is set to 
zero because the bound state is slightly off-shell and free from IR divergence. The brackets 
(■ ■ ■) mean that we take the difference between the expectation values with respect to the 
spin J = 1 state and the spin J = state. R' is the two or more photon exchange part R 
of Green's function Gq given by Eq. ([TT| ) minus the pole contribution of the ground state: 



^'-.-J^-^)- (-' 
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The explicit form of R' can be found in Eq. (28) of Ref. [Q. Integration over / and g in (|A6| ) 
can be carried out using the identity for R' given in Ref. 



where 

i?(g)=ln2-| + ^tan-i(i)-iln(l + 4)+^^. (A9) 

2 q V7/ 2 V 7^/ g^ + 7^ 

Unhke the two or more photon exchange part, the zero and one photon exchange parts 
of the Green function give rise to hnear and logarithmic divergences, respectively. These 
divergences must be taken care of by the NRQED contact terms. Thus it is convenient to 
treat the calculation of the bound state and the corresponding contact terms together, since 
they have the same UV behavior. 

For the one photon exchange part the bound state calculation gives 

AE[Fig.2(b)] = 2 / ^ ',, ,,-,, :1'L ■,.^. (^Hr,0) 



(27r)12 (p2 + ^2)2(^2 + ^2)2 ^ - V ' ^' f2 + ^2 

2m 

g^ + f 



xVcil,q)^^—-^VK{q,p) 



-21n2 + 21n(-')-^ 
7/ 16 



(AlO) 



where A is the UV cut-off. It is put to infinity for the finite parts. 

The corresponding contact term should have the same structure and we find it among 
the two-loop contact terms. It is one Coulomb photon exchange between the Fermi and the 
relativistic kinetic energy potentials. (See Fig. ^ (d).) Since this is a contact potential, the 
first order perturbation theory with the potential becomes the wave function at the origin 
squared |0(O)p = '-f^/ir times the minus sign of its scattering amplitude. To avoid the IR 
singularity caused by the on-shell condition, we put the infinitesimal photon mass A in the 
photon propagators. This leads to 



AE[Fig.2(d)] = |0(O)|2(V2e-k(i)) = Ep(Z«)2 [-21n(^) + 21n3 

where V2c-k{i) is given by (^61) . This In A cancels that of Aii^[Fig.|^(b)]. 
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(All) 



The zero-photon exchange part of the bound state theory can be calculated in a similar 
way. (See Fig.0(c).) One subtle thing takes place in the calculation of the zero photon 
exchange part, because linear divergence causes the integrated result to depend on how it 
is parameterized. If the regulated momentum is shifted, it may give an additional finite 
piece to the answer. In order to get rid of the uncertainty due to linear divergence, we will 
consider the bound state calculation and the contact term calculation together and put the 



cut-off A in the corresponding transverse photon line p2| . We choose the momentum k of 
the transverse exchanged photon line to have the cut-off A throughout the whole calculation, 
since it is easily recognized. Keeping the cut-off A at finite value, we get 



d^p f^ d^k 



A^Ps-2Wl = 2/75^/ 



iVrry 



^^2 



X VFip + k,p) ^^ 2 ^x(p) 



EpiZa) 



The corresponding contact term gives 



7' 



77r 8 



(A12) 



AE[Fig.2(e)] = ^(Vic-k) = EF(Za) 



n 



A' 
-2 — 

77r 



(A13) 



where Vic-k is given by (|25|) . In these calculations, we have used the expansion 



tan 



if^ 



71 X 1 /A 



+ 



X 



(A14) 



XJ 2 A 3VAy ' VVA. 

Clearly, the one-loop contact term kills the linear A divergence in A£'[Fig.0(c)]. 

The other two-loop contact terms are linearly divergent due to the threshold singularity. 
The diagram Fig. ^(f ) gives 

.'3 r A _ 2 



7 



AE[Fig.2(f)] = ^(V2c-k(2)) = EF(Za 



n 



-4— -21n- + 3 

ATT 6 



(A15) 



where V2c-k(2) is given by ([A4|) . The diagram Fig.H(g) has an expression different from the 
diagram Fig. |](f), but their integrals are identical: 

-'3 r A _ 2 



AE[Fig.2(g)] = ^(V2c-k(3)) = EF(Za)^ 



-4- 



Att 



2 In - + 3 



(A16) 
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where V2c-fe(3) is given by ([XSl) . 

The remainder is the contact term that consists of a one-loop contact term and one 
Coulomb photon exchange: 

.3 r A 



A^[Fig.2(h)] = ^(V2c-k,iioop2) = EF(Za) 

TT 



Att 



where V2c-k,iioop^ is given by 

Summing up all the contributions relating to the kinetic term correction, we get 



(A17) 



AEk = EpiZaf 



i_61n- + 21n('- 
2 3 V7 



(A18) 



Each linearly divergent diagram could have a value which depends on the regularization 
method, but the sum remains the same because of gauge invariance. 

3. The Darwin Term 

The computation method of the Darwin term shown in Fig. ^is almost identical with that 
of the kinetic term. The bound state calculation is also from the second order perturbation 
theory. The number of graphs has been already taken into account in the following equations: 



Ai?[Fig.3(a)] = 2| 



d^pd^qdHd? 



{8V^r 



(27r)12 (^-2+ ^2)2(^2 +^2)2 

^ /-647r\ ~ - 



= EpiZa) 
Similarly, we find 

AE[Fig.3(b)] = EpiZaf 
AE[Fig.|(c)] = EpiZaf 



ZaY 
3 
'2 



-\n(-] +ln2 

A 

+ 

'-fir 



(A19) 



(A20) 
(A21) 



The contact terms give the following contributions: 
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AE[Fig.3(d)] = -2^ / -^(l^^(0,A:))^^\/z.(A:,0) 



EpiZaf — 

77r 



A^[Fig.|(e)] = EF{Zaf 
A^[Fig.|(f)] = EF{Zaf 
AE[Fig.|(g)] = EpiZaf 



A 



A^[Figi(h)] = Ef{Z 



a] 



A 2 

2— + ln--l 
. Avr 3 

^ A 2 3 

2—- + In 

Avr 3 2 

A 

-4— + 3 

ATI 



The sum of the contributions ( [A19|) - (|A26|) is 
A£;, = EpiZaf 



■„^U3,4 



(A22) 
(A23) 
(A24) 
(A25) 
(A26) 



(A27) 



4. The Derivative Fermi Term 



The corrections involving the derivative Fermi term shown in Fig. ^ are calculated in 
the first order perturbation theory. We obtain 



A£;[Fig.4(a)] 



(27r)3 7 (27r)3(|p'_A;|2+y)2(fc2+^2)2 

EpiZa)^ 



{Vw{p,p-k)) 



-2A + ^ 
77r 2 



AE[Fig.|(b)] = EFiZaf2 — 

77r 



AE[Fig.|(c)] = EpiZaf 



AE[Fig.^d)] = Ef{Z. 



a] 



AE[Fig.|(e)] = EpiZa) 



The sum of the contributions 



4A + 21n^-3 
Avr 3 

A 2 

4— + 21n--3 
Avr 3 

A 

-8— + 6 

Att 



-([M^is 
AE^ = EF{Zaf 



.,2 5 

4 In - + - 

3 2 



(A28) 
(A29) 
(A30) 
(A31) 
(A32) 



(A33) 



30 



5. The Seagull Term 



The diagrams involving the seagull term are shown in Fig. ^ The bound state calculation 
is also carried out in the first order perturbation theory. The results are 

d^p /-A d^k f dH (8x/7rp)2 



AE[Fig.5(a)] = 2 / 



(27r)3 7 (27r)3 7 (27r)3 (|p'_ A;|2 + 72)2(1^ _ /|2 + y)2 
X {Vsip-tp-kJ)) 
= Ep^Zaf 



AE[Fig.|(b)] = Ep^Za) 

AE[Fig.|(c)] = EpiZa) 
AE[Fig.|(d)] = Epi^Za) 



A_inf^)+ln2 

77r \7 

_ A" 

77r 

-4A + in('^')+| + ln3-21n2 
Att V a / 2 

4A-3 

Att 



(A34) 
(A35) 

(A36) 
(A37) 



The sum of the contributions ( |A34|) - (|A37|) is 

'A^ 



/\Es = EpiZa) 



-In 



In 



2 1 



.7/ 3 2 

Adding up Eqs. ( |AT8| ), (|A27D , ( [A3^ ), and (|A38D we obtain the Breit {Zaf term 



(A38) 



A^Breit = 2^^(^")'- 



(A39) 



APPENDIX B: CALCULATION OF THE a{Zaf CORRECTION 



1. Vacuum Polarization Contribution 



First let us calculate the NRQED corrections related to the vacuum polarization inser- 
tion in the transverse exchanged photon. The bound state calculation is the first order 
perturbation theory with the NRQED potential Vtvp of (|38|): 

Tc a f d^p [^ d^k 



A^[Fig.6(a)] = Ef- 
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a{Zay 8 



2 m^ 



n 15 



A _ 3 

77r 2 
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Similarly, for the scattering diagrams of Fig.|](b) - (e), we obtain 

a(Z«)2-8 A 



AE[Fig.6(b)] = Ep- 
AE[Fig.|(c)] = Ef 
AE[Fig.|(d)] = Ef 
AE[Fig.|(e)] = Ef 



n 15 77r 



a{Zay 8 
n 15 

a{Zaf 8 
n 15 
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^ 15 
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2 17 

-2- 21n- + — 

Att 3 12 

A 3 17 

-2 31n2 + -ln3 + — 

Att 2 6 

A 15 
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Next let us calculate the contribution coming from the NRQED potential representing 
the vacuum polarization insertion in the Coulomb photon given by the potential Kvp of Eq. 
(PSp. The bound state calculation is in the second order perturbation theory. The two- or 
more-photon-exchange part of the Coulomb Green's function gives the correction 



AE[Fig.6(f)] = 2Ef 
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Similarly, the one-photon-exchange part gives 

a{Zay 8 



AE[Fig.6(g)] = Ef- 



The zero-photon-exchange part gives 
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71 15 
The corresponding scattering diagrams give the following corrections: 

-4:7r{Za)a f^ d^k k^ -2m k^ 



AE[Fig.6(i)] = -2Ef 
= Ef 
Similarly, we obtain 
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A£;[Fig.6(j)] = Ef 

AE[FigJ|(k)] = Ef 

AE[Fig.||(l)] = Ef 

A^[Fig.|(m)] = Ef 
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The sum of AE[Fig.|(a)] ■ ■ ■ AE[Fig.|(m)] gives the NRQED contribution of AEnrqedvp 
ofEq. §3). 

The remainder is the QED scattering diagrams contribution. The QED three-photon- 
exchange skeleton diagram gives the hyperfine splitting contribution: 

(Za^ f°° /•°° /-i lQ{p^ -p- q + q'^) 
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/•oo r I 
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Thus the vacuum polarization insertion in the middle exchanged photon gives 
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AE[Fig.6(n)] = Ef^^ / dtp2{t) / dp dq d{cost 
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where p2it) is the second-order photon spectral function given by 
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The vacuum polarization insertion in the outermost exchanged photon leads to 
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16(p2 —p. q-\- g2^ 
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The scattering amplitude due to the KP contact term is 
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(p2 + A2)(p2 + 4^2(1 _ t2)-l)(^2 + A2) ' 

The sum of AE[Fig.§(n)], AE[Fig.§(o)],and AE[Fig.|(p)] gives the QED contribution 
A^QEDvp of (ill). 
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2. Radiative Photon Contribution 



We list formulas needed for the calculation of the bound state contribution Eb- For the 
diagram Bqq of Fig. |^ we find 



AE[Fig.7(Boo)] =2Ef— — 
= Ef 



TT 2a f d^pd 



ly/rn 



^^2 



P-P 



dk- 



-2m 



7^ StT J (27r)6 (p2 _|_ ^2~)2(^^2 _|_ ^2~)2 ^ J^ p2 _^ p^2 
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3 7 7o (po + 7)^ 

Some contributions of Fig. |^ can be similarly reduced to integrals over the radiative photon 
momentum k while others are harder to simplify: 
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AE[Fig.KBn)] = 2Ef 
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For the scattering state, we put 
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The NRQED scattering state contributions to Es are given by the following terms: 
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Similarly, we find 
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The g^ derivative Fermi term A£^g2 contributing to Eji is calculated by replacing ff +p''^ 
in VV of (^2|), etc. by g^. We obtain 
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AE[Fig.|(c)] = Ep{Zaf 
AE[Fig.|(d)] = EpiZaf 
AE[Fig.|(e)] = EpiZa)^ 
The sum of the contributions (|B39|) - (p343|) is 
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(B44) 
we obtain the contribution 



Multiplying AEq2 with the "renormalization" constant Cg2 of ( 
ASq2 listed in Table 0. 

The reminder is the contribution ASp/p coming from the p'p coupling interaction. The 
p'p potential Vpip{p',p) for spherical symmetric states is defined by 
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Vp'p{p\p) 



Zan —p'^q'^ + (p ■ q)"^ 



With this expression, the bound state calculation of Fig. ^(a) gives 
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The 2-loop scattering diagram of Fig. ^(b) gives the contribution 
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The sum of (|B46|) and 
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where A is the infrared cutoff. Multiplying this with the " renormalization" constant Cp/, 



of 



we find the contribution /S.8pip listed in Table III. 



The contribution Eq of the QED scattering diagrams shown in Fig. |T^ can be calculated 
using the techniques similar to one described in Ref. 0. These diagrams may be expressed 
in the form: 
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where £fj,u\ and M.^^'^^ are factors representing the electron line and the muon line, respec- 
tively. For instance, £^u\ corresponding to the diagram Tia of Fig.|^ is given by 



{{p + ly -w? + ^e)((? + ly -m? + i(^){{k + q + ly - m'^ + ie){{k + ly - w? + ie) 
where / = {m, 0). M.^^^ may be written as the sum of six permutation terms: 
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where r = (M, 0). 
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The integral is greatly simplified in the limit of infinite muon mass. We can extract 
the contribution to the hyperfine splitting from each diagram using the projection operator 
given by Eq. (65) of 0]. For instance, the contribution from the diagram Tia of Fig.|10|, after 
carrying out the k integration and subtracting the vertex renormalization term, is expressed 
with the help of Feynman parameters Zi, Z2, and z^ as 

AE[Fig.lO(Tia)] = 2Ef^ ^-- / dz,dz2dz,5{l - zi - z^ - z,) 

d^p / d^q- 
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where 



^1 — Z5, Aiq — 1 — Zi, A2q — ~Zi, 

V = Zi+ Z2- {Zi + Z2)Ai + Zi Aiqq^, 

Vo = Zi + Z2- (zi + Z2)Ai. (B53) 

The integral A£'[Fig.|^(Tia)] has one threshold singularity at g = and another at p = 0. 
The threshold singularity at g = is canceled by that of the NRQED scattering diagram 
Aii^[Fig.[iy(A)] which consists of one Fermi and two Coulomb potentials multiplied by the 
Fermi-term " renormalization" constant of the Fermi term Og, namely the anomalous mag- 
netic moment of the electron. The latter is of the form 

a{Zaf 1 



AE[Fig.lO(A)] = -2Ef^^^^-- j dz^dz2dz^5{l -zi~Z2- z^) 



-I^'pI^'o r-m^ . A2WU- ' -1 . , yu.. , ,./ -(->^Al + 8^4 . (B54) 
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The threshold singularity at p = is canceled by the contact term consisting of the lower 
order contact term (the KP term), which contributes to the a{Za)EF correction, and the 
Coulomb potential. /,From the second and fourth diagrams of the KP diagram of Fig. 
p!0|(KP) we find that the contact term contribution corresponding to Tia and A of Fig. |10| is 
given by 

1 
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71 071^ J 
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The sum of ( [B52D , ( |B54| ) , and ( [B55D still suffers from a severe infrared singularity in the limit 
of vanishing radiative photon mass. In order to perform numerical integration we identified 
the IR singular terms of (|B52| ) and ( [85 51) and subtracted them lirora each integral. The IR 



subtraction term for (P52| ) is of the form: 

a(Za^2 
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where 
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IR 



zi + Z2- {zi + Z2)Ai + ziq"^ . 



(B57) 



This IR subtraction term is completle canceled by that for Tq of Fig. ^. Similar cancellation 
occurs among the diagrams Tn,, T2, and T3. The IR subtraction term for Eq. ( P55|) is given 
by 
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This type of IR sungularities of the KP contact terms cancel out completely among them- 
selves. When summed over all diagrams of Fig. ^, the resultant integrand has only the 
infrared singular terms of the form ln^(A/m) and ln(A/'m). 
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FIG. 1. Comparison of the QED and NRQED scattering amplitudes. The small circles on the 
edges of each diagram indicate that a diagram is evaluated on-the-mass-shell and at the threshold 
external fermions. The shaded circle in (a) is a contact term. The doubly-shaded circle in (b) is 
a contact term, too. A wavy line is covariant photon in the Feynman gauge, a early and dashed 
line are transvers and Coulomb photons, respectively, in the Coulomb gauge. The muon is treated 
as the static external source and indicated by a dot. For NRQED electron vertices, consult the 
NRQED Feynman rules given in Ref. [3] . 
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FIG. 2. Relativistic kinetic term diagrams. The brackets ( ) indicate that the diagram is 
evaluated with the bound state wave function. The contact term diagram is shown in the right 
below of the corresponding bound state. The time-reversed diagrams are not shown. R in (a) 
represents two or more photon exchange part of the Coulomb Green's function. 
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FIG. 3. Darwin term diagrams. 
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FIG. 4. Derivative Fermi term diagrams. 
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FIG. 5. Seagull term diagrams. 
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FIG. 6. Vacuum polarization diagrams contributing to the a{Za)'^ correction. 
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FIG. 7. Diagrams relevant to the NRQED bound state contibution Eb of the a{Za)'^ correc- 
tion. 
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Diagrams relevant to the NRQED scattering state contibution Es of the a{Za) 
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FIG. 9. Diagrams relevant to the "renormalization" constant contibution Eji of the a{Za)'^ 
correction. Only the diagrams involving p'p interaction term are shown here. For other terms, see 
Figs. 2,3,4, and 5. 
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FIG. 10. Diagrams relevant to the QED scattering state contibution Eq of the a(Za)^ correc- 
tion. The last line is the difinition of the KP potential which contribute to the a{Za) correction 
of the hyperfine splitting. 
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TABLES 
TABLE L Bound state contibutions Eb of the a{Za)'^ correction from individual diagrams 
shown in Fig. 7. The overall factor Epa(Za)'^ /tt is omitted. 



Diagram 



Contribution 



5, 



00 



16 
3 



2Am 

^,2 



fin 



2Am 



+ fln2-| 



5, 



01 



16 
3 



1 1^2 /^ 2Am 
4^11 I 2 



ln2-i)ln(^| +iln2 + ln=^2 + 4 



fi, 



02 



16 
3 



A _ 49 _ 



i3in(2^j_||_^ln2 



B 



10 



iln=^ ^ 



31n2 + i)ln(^) +In2 + 31n22 + ^ 



Bu 



81 ll 
3 I 6 



In 



2Am 

7^ 



1-10.499 478 (61) 



B 



12 



3 \ 6 ~'~ 4 



In 



2Am 



- 1 - 5.393 042 (55) 



B^ 



20 



+0.761 648 (52) 



B^ 



21 



-0.491 192 (42) 



B. 



22 



+0.221 861 (36) 



B 



do 



ilnfMlli] _ iin2- A 



24 



B. 



dl 



-0.316 379 (13) 



B. 



d2 



-0.116 135 (19) 
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TABLE II. Scattering state contributions Eg of the a{Za)'^ correction from individual diagrams 
shown in Fig. 8. The overah factor Epa{Za)'^ /it is omitted. -Lz(l/3) = 0.366 213 229 • • • is the 
value of dilogarithmic function Li{x) at x = 1/3, where Li{x) = — Jq dt '^^ ' . 
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TABLE III. Renormalization constant contributions Eji of the a(Za)^ correction from indi- 
vidual diagrams shown in Fig. 9. The overall factor Epa{Za)'^ /it is omitted. 



Diagram Contribution 
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31n|-lnf^Y 

3 \lj 
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TABLE IV. The QED scattering state contribution Eq for various photon mass A. The 
corresponding diagrams are shown in Fig. 10. The overall factor Epa{Za)'^ /tt is omitted. The 
uncertainty of Eq is from numerical integration by VEGAS. 



AVm2 


Eq 


Eq - {EQ)\og 


uncertainty in Eq 


10-5 


-90.075 1 


-8.697 1 


0.0024 


10-5.5 


-107.943 8 


-8.707 8 


0.0033 


10-6 


-127.575 1 


-8.713 9 


0.0040 


10-6.5 


-148.971 3 


-8.717 6 


0.0048 


10-7 


-172.139 


-8.725 5 


0.0029 


10-7.5 


-197.062 7 


-8.722 


0.0030 


10-8 


-223.758 9 


-8.723 7 


0.0026 


10-9 


-282.449 4 


-8.723 4 


0.0016 
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